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1 Introduction 

A remarkable characteristic of the general-relativistic theoretical framework, encompassing 
classical as well as semiclassical aspects, is that it seems to contain the seeds of its own 
destruction. We shall argue that, when properly interpreted, this framework also incorpo- 
rates clues for efficiently debugging its problems. For instance, classical general relativity 
predicts the unavoidable formation of black holes through gravitational collapse. But then, 
the classical laws of black hole dynamics call openly for quantum mechanics for consistency. 
A crucial step towards this consistency was Hawking's derivation that black holes should 
emit thermal radiation, and therefore evaporate, due to quantum effects [1, 2]. However, 
right after Hawking presented his derivation, Unruh showed that it was somewhat un- 
trustable due to what has come to be known as the trans-Planckian problem [3]. Hawking 
radiation appears to come from the red-shifting by the collapsing geometry of modes with 
huge initial frequencies, well beyond the Planck frequency ujp. To estimate how huge these 
frequencies are, one only has to calculate the blue-shift factor e KH * BH , where kh stands 
for the surface gravity of the initial black hole («h = (4m) _1 , with m its mass) and tpH 
is the lifetime of the evaporating black hole (ipH = 51207r(m/mp) 3 £p, where mp and tp 
are the Planck mass and time respectively), and multiply it by the typical frequencies of 
Hawking radiation, of the order of kh (for the argument that follows, it is not necessary 
to take into account the variation of kh during the evaporation). For a Solar mass black 
hole (kh — 10 _39 wp and ipH ~ 10 118 tp) this calculation yields a value 

^max ~ K H e KHtBH ~ 10 lo79 o;p (1.1) 

for the highest invoked frequency w m ax- Just one second after the ignition of the black 
hole [4] the highest frequencies involved in the Hawking process are already of the order 
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w max ~ 10 up. It is extremely difficult to believe in the validity of quantum field 
theory in a curved classical background up to these enormous energy scales. As Unruh 
used to say, this should be considered the trans-Planckian disaster or, perhaps, the trans- 
Planckian catastrophe, in analogy with the ultraviolet catastrophe of the late 19th century. 
In other words, the evaporating black hole scenario based on this mechanism for Hawking 
radiation can be considered strongly unsatisfactory. 

At these extremely large energy scales the geometry itself should be "quantum" (using 
the word in a vague sense) and modify the effective behaviour of the matter fields. If one 
could prove that Hawking emission is not really dependent on these effects of the unknown 
physics at huge energies, then the evaporating black hole scenario could turn out to be 
reliable after all. 

The assumption that long-lived trapping horizons exist in the low-energy physical 
description (a prediction of standard general relativity) calls for some hypothesis about 
the characteristics of the high-energy physics. It is difficult, if not impossible, to catalogue 
all the possible scenarios that might be at work at high energy scales. However, it seems 
reasonable that they could be separated into two classes: those in which the classical notion 
of a long-lasting trapping horizon is preserved and those in which it is not. Does the causal 
unidirectional property of a trapping horizon survive to high-energy physics corrections or 
is it only a low-energy property? The simplest way to grasp the effects associated with 
these two possible scenarios is to analyze linear effective quantum field theories containing 
dispersive corrections of either subluminal or superluminal character, respectively. We 
cannot logically discard the possibility that the underlying high-energy theories, which 
either preserve the horizon completely or maintain it only as a low-energy property, do not 
give rise to the predictions of these sub/super luminal effective field theories. However, we 
find this possibility difficult to imagine and, in principle, unjustified, as it would contradict 
the very essence of using effective physical descriptions. 

These dispersive analyses were introduced by Jacobson in [5, 6] and soon followed 
by Unruh [7]. At present there are many works analyzing the effects of non-relativistic 
dispersion relations (see the review [8] and references therein). Typically the conclusion 
of these investigations can be summarized by saying that Hawking-like emission survives, 
under mild conditions, the presence of modified dispersion relations of both kinds, and so 
that the evaporating black-hole scenario is reliable. Sometimes, it is even tacitly assumed 
that the real (astrophysical) Hawking radiation would be produced in one of these modified 
manners. However, we want to highlight in this paper that one cannot just study the effect 
of the modified dispersion relations in Hawking emission without further considerations. 
Rather, one has to probe the consistency of the whole scenario as well, analyzing collateral 
effects of dispersion, which may well become central. 

In this paper we shall present a comparative discussion of the three different scenar- 
ios that one might envisage in the light of the trans-Planckian disaster of the standard 
evaporating black hole scenario. The stress will be put on how a final consistent scenario 
might look like. In the first two scenarios it is assumed that a long-lived trapping hori- 
zon has been formed from the point of view of the low-energy physical description. The 
first scenario (Section 3) is one with subluminal dispersion relations: the notion of horizon 
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remains in the underlying high-energy theory. The second scenario (Section 4) involves 
superluminal dispersion relations: the notion of horizon does not survive in the underlying 
high-energy theory. Finally, in the third scenario (Section 5) we assume from the start that 
long-lived trapping horizons do not form in the first place. Black holes are substituted by 
black stars, bodies hovering right outside but extremely close to their gravitational radius. 
We will show that there exists a mechanism by which these bodies can emit Hawking-like 
radiation. In this mechanism huge trans-Planckian frequencies play no role whatsoever, so 
it can be perfectly described with unmodified relativistic dispersion relations. 

2 Preliminaries 

The original Hawking derivation of the thermal emission relies strongly on the dynamical 
formation of a horizon. It is similar to a process of "cosmological particle production" 
but modified due to the presence of a horizon. Cosmological particle production is tied up 
to the dynamics of spacetime: when spacetime stops expanding and becomes stationary, 
the production also stops. However, in the Hawking process, even though the dynamics 
is the agent responsible for the particle production, the particle production never ends 
once a horizon has formed (in the absence of back- reaction) , even if the dynamics stops. 
The infinite red-shifting effect of the horizon makes the short movie of events right before 
the formation of the horizon to be contemplated from the asymptotic region as a never- 
ending movie. In the asymptotic form of the radiation only the surface gravity of the black 
hole enters into play, controlling the temperature of the thermal spectrum. However, as 
explained in the introduction, this huge red-shifting effect is in turn responsible for the 
trans-Planckian problem. 

With modified dispersion relations starting to show up not far from the Planck energy 
scale, the dynamics of the formation of the horizon turns out to be unable to produce 
a stationary emission of particles from the horizon. This is so because subluminal and 
superluminal behaviours of the types described below produce an effective cut-off in the 
frequencies. But, as is well known, a plain Planckian cut-off Co> max completely kills Hawk- 
ing radiation [5, 9]. However, the presence of modified dispersion relations only kills the 
standard dynamical way in which Hawking radiation appears, but opens up the door to a 
quite different mechanism of particle production which does not rely on the dynamics but 
rather on the characteristics of the configuration itself. This new mechanism, known as 
mode mixing, occurs whenever, for a particular positive energy (a conserved quantity in 
stationary configurations), there exist modes with both positive and negative norms with 
respect to the field theory inner product (see for instance 2.2) below). 

An elegant way to treat the appearance of dispersion is the covariant formalism devel- 
oped by Jacobson [10]. Apart from the metric field, there is a unit timelike vector field u a 
singularizing the specific reference frame in which one has to measure the scale at which 
the modifications appear. Using for our discussion the simplest case of a scalar field, the 
field theories we are interested in satisfy the equation 



□ + F(h ab V a V b ) = 0, 




- 3 - 



where h = g + u a u . In the space of solutions of the previous equation one can define 
a conserved inner product as 

(01, 02) = ~i J d&z a (<t>\ d a fa) . (2.2) 

In the reference frame defined by the vector field u a the dispersion relation associated with 
this equation in the slowly-varying approximation reads 

J 2 = k 2 + F{k 2 ), (2.3) 

where k 2 = h a \ ) k a k h . 

3 Subluminal dispersion 

We will consider two types of subluminal dispersion. In the first one, to 2 is monotonically 
growing with k and approaches a constant value for k — > oo. In the second one, uj 2 increases 
from k = till it reaches a maximum value and then decreases back to zero at a finite k. 
Two specific functions representative of both types of behaviour are 

F{k 2 ) = -k 2 + fc|tanh(A; 2 /fcp) , F(k 2 ) = -k A /kl , (3.1) 

where kp is the Planck wavenumber. In particular, throughout our discussion we will for 
simplicity assume spherical symmetry and use an effectively one-dimensional spacetime 
with a metric of the form 

ds 2 = -[1 - v(x) 2 ]dt 2 + 2v(x)dtdx + dx 2 , (3.2) 

where v is the free-fall velocity from infinity, i.e. the flow velocity in the language of acoustic 
metrics (we have set the speed of light c = 1). In addition we shall choose u a = {l,f}- 

As shown initially by Unruh [7], mode mixing in the surroundings of the horizon 
produces a stationary spectrum of particles with thermal properties, if uj c 3> kh, where ui c 
is the threshold frequency below which mode mixing takes place. The spectrum is perfectly 
Planckian up to frequencies of the order of lo c . It should be noted that some subluminal 
models still exhibit a trans-Planckian problem. The first subluminal dispersion in (3.1) 
with asymptotically vanishing flow velocity v (x — > oo) — > is such an example: finite-A; 
Hawking modes originate via mode conversion outside the horizon from infinite-/c modes 
in the past, see e.g. [11]. 

Therefore, in general: yes, a subluminal mechanism could be responsible for black hole 
evaporation while avoiding the trans-Planckian problem. However, this scenario as a whole 
maintains some other problems of the standard evaporation picture. On the one hand, it 
is well known that the existence of trapped surfaces entails the formation of singularities 
in standard general relativity. Any new theory beyond standard general relativity which 
incorporates a subluminal behaviour would maintain this tendency to form singularities. 
Thus, in this scenario one would have to face the information- loss problem [12]. It is 
difficult to see how all the information that went to form a black hole could be released in 
the very last stages of its evaporation. 
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On the other hand, the evaporation of a stellar mass black hole would take 10 times 
the current age of the universe. The existence in this scenario of long-lived trapping 
horizons enormously delays the time at which we will be able to have experimental feedback 
about the real nature of the underlying high-energy theory. Therefore, in the absence of 
primordial black holes [13] or microscopic black holes at the LHC [14] 1 , this scenario does 
not allow an open positive exploration of the full nature of its underlying physics, at least 
not to the human race as we understand it. As physicists we find this idea, to say the least, 
disturbing. 

4 Superluminal dispersion 

Superluminal scenarios are completely different. The physics at work inside the horizon 
becomes exposed to the outside world. In particular, the characteristics of their internal 
geometry can have an influence on how black holes radiate and more in general on how 
they behave. The analysis of stationary configurations with a black hole horizon has 
shown that mode mixing is also able to provide a mechanism for Hawking-like emission in 
a superluminal scenario (see [15] and [16] for a more recent approach). In those analyses, 
there is always a crucial assumption regarding how the in-vacuum state looks like near the 
place where the internal classical singularity is supposed to be. This amounts to selecting 
a boundary condition at the singularity. One way to have a good control of this boundary 
condition is to substitute the place where the classical singularity would appear by an 
internal asymptotic region. Different vacuum states can be selected by allowing or not the 
net entrance of particles from this internal asymptotic region. 

Contrarily to subluminal scenarios, any new theory beyond standard general relativity 
that incorporates a superluminal behaviour would blur the very essence of the horizon. 
There will be high-wavenumber signals that can escape from the hole. Actually, assuming 
that the superluminal correction does not saturate, at any location inside the black hole, 
there will always be a wavenumber above which the signals can escape towards the exterior. 
Under these circumstances it is no longer sensible to assume that the theory will have a 
tendency to form singularities. Quite the opposite view would be more reasonable, namely, 
that the internal region is perfectly regular so that any wave entering into the hole would 
after some time travel back towards the horizon. Superluminal dispersion would therefore 
not only remove the horizon, but also smooth away the curvature singularity. 

When in the internal region one allows the existence of particles travelling towards 
the horizon, that is, when the vacuum state is not the so-called in vacuum, the particle 
spectrum in the asymptotic region could be largely distorted from a Planckian form (see [17] 
for a discussion on the various vacuum states one could set the system in). For instance, 
a similar situation can be seen in the analysis of configurations with a black and a white 
horizon where a black-hole-laser effect is generated [18-21]. There, the existence of an 

1 We find this possibility much more speculative. In turn, the existence of these microblack holes would 
be revolutionary, implying that the relevant quantum gravity scale is not the Planck scale but the much 
smaller electroweak scale. 
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internal white horizon completely destroys Hawking radiation, giving place instead to a 
self- amplified creation of particles. 

Moreover, as we are going to show, under reflecting boundary conditions at some point 
located inside the hole, the theory becomes unstable, as already observed in the context 
of BECs [22, 23]. This means that there appear frequencies in its spectrum with positive 
imaginary part (again, very much as in the presence of a combination of a black-hole and 
a white- hole horizon). To be more specific, let us analyze a scalar quantum field with a 
quartic superluminal dispersion relation 

F{k 2 ) = k 4 /4 (4.1) 

in a one-dimensional spacetime with a metric (3.2) in which the velocity v has a step-like 
profile, i.e. such that it vanishes in the exterior and has a constant and homogeneous value 
v > 1 in the interior. As before, the spatial dimension represents the radial coordinate of a 
more realistic spherically symmetric three-dimensional configuration. Furthermore, let us 
assume that the field is fully reflected at some point inside the hole located at a distance L 
from the horizon (corresponding to an s-wave analysis and a symmetry condition at r = 
in three dimensions). Then, from the numerical analysis, we can estimate the number of 
instabilities nj ns t, as well as their maximal typical growth rates r max and the real part oj n 
of their frequencies, as functions of v and the size L of the inner region. The details of 
this numerical instability analysis can be found in the appendix. Overall, the qualitative 
behaviour of these quantities is 

r max ~ 2(v - 1)/L, n inst ~ 2w c /r max , (4.2) 
uj n ~ ±r max n (n = 1 . . . n inst /2), (4.3) 

where uj c is the critical frequency (frequency above which there is no more mode mix- 
ing [16]). These results are in agreement with the numerical and analytical studies carried 
out for the black/white-hole case in ref. [19]. The interpretation of these quantities is 
straightforward and could actually have been obtained from a qualitative analysis of the 
problem, as discussed in what follows. r max is directly related to the inverse of the time 
needed to travel back and forth inside the hole. The maximum real part of the frequency of 
these instabilities is just oj c . These two observations are closely related to the fact that the 
leaking of resonant modes inside the hole is responsible for the instability. The number of 
instabilities is determined by the ratio between these two quantities w c /T max , at least for 
sufficiently large L. Again, this is just a consequence of the fact that the resonant modes, 
which obey a kind of quantization rule, build up and eventually leak through the slightly 
permeable horizon giving rise to the unstable behaviour. 

As an example, figure 1 shows the map of instabilities for a black hole whose interior 
size is 40fcp 1 and the flow velocity is twice the speed of light. In this plot, the growth rate 
Im(cj) is measured in units of r max and Re(o;) is measured in units of u c . 

If we consider a Solar mass Schwarzschild black hole written in Painleve-Gullstrand 
coordinates (in which the radial part of the metric is of the form (3.2)), it is a trivial 
exercise to show that the average flow velocity in the interior is precisely twice the speed 



- 6 - 




Figure 1. Map of instabilities for a black hole with superluminal dispersion and reflecting boundary 
conditions, whose interior size is ^tOfcp 1 and with a flow velocity twice the speed of light. The growth 
rate Im(u) is measured in units of r max and Re (a;) in units of ui c . 



of light. Since its size is typically 3 km, we have the following numbers for the quantities 
above: 

r max ~ 2 • 10V\ n inst ~ 10 38 , w c ~ 10 43 s- 1 . (4.4) 

The essential ingredients for this result are the appearance of additional propagation chan- 
nels characterized by a high-energy scale (kp in (4.1)) and a superluminal character. Our 
qualitative conclusions are therefore valid beyond the particular model that we have an- 
alyzed numerically. Note that although the temporal scale of the instability is very large 
compared with Planck time and would therefore be seen as metastable from the artificial 
black-hole perspective, it is actually tiny in astrophysical terms. This means that a Solar 
mass black hole with superluminal dispersion would be extremely unstable in astrophysical 
terms and would tend to eliminate its horizon in microseconds. Thus, we see that if the 
true curvature singularity is smoothed away due to the presence of superluminal dispersion 
relations (as discussed above), black holes as such would be very unstable. In other words, 
this analysis suggests that any theory beyond classical general relativity that incorporates 
superluminal behaviour will avoid the formation of long-lasting trapping horizons. 

5 Avoidance of long-lasting horizons 

The last scenario that we will consider is one in which we assume from the start that 
there are no astrophysical bodies with strict long-lasting trapping horizons. Black holes 
are substituted by material bodies with a real surface which is hovering right outside, 
but extremely close, to their gravitational radius. The existence of these objects seems 
largely implausible within standard general relativity, but they might occur in a modified 
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theory of gravity such as semiclassical general relativity (see the discussion in [24]) or 
condensed-matter-inspired emergent gravity [25]. Using the terminology of [24] we will 
generically call them black stars. There are several black-hole mimickers proposed in the 
literature [24, 26, 27]. What we are going to describe here does not depend on the specific 
reasons behind their existence or on their internal constitution, but only on the fact of 
staying close to horizon formation. 

The entropy associated with these black stars would not have a purely geometrical 
origin as is typically assumed for black holes. Instead, it will be associated with their matter 
content under the influence of gravity. It is however remarkable that there exist strong 
indications that the entropy associated with matter configurations on the verge of forming 
a horizon would follow an area law [28-32]. It is also remarkable that objects of this kind are 
able to emit radiation by a process closely related to the one originally devised by Hawking. 
In [33] it was shown that objects without horizons can emit Hawking- like radiation during 
any period of exponential approach to horizon formation. This observation is resonant 
with the fact noticed previously in [34] that a collapsing shell bouncing or stopping close 
to the horizon will emit a burst of thermal radiation. 

For instance, imagine a black star vibrating between two radii rn + x ut and rn + Xm, 
with both x ou t/rH, Xi n /rH <C 1. To characterize the particle emission associated with a 
vibration cycle we can use the techniques developed in [4, 35, 36]. Let us consider a 
simple scalar quantum field over a 1+1 geometry. The amount of particle production 
will be characterized by the Bogoliubov coefficients, which can be calculated from the 
expression [33] 

/W = I due-^VMe-*™, (5.1) 

provided that the relation between two different outgoing null coordinates of the geom- 
etry U = U(u) is known. Typically, U and u represent asymptotic past and future null 
coordinates, respectively. Instead of using a complete model for the vibrating geometry, 
one can alternatively use a stationary Schwarzschild-black-hole geometry and encode the 
vibrational features into the definition of a non-stationary vacuum state (see [4] for a more 
detailed discussion). With this aim, let us take a vibrating timelike trajectory r(r) con- 
necting t"h + x ou t and rn + Xm- An asymptotic observer will label the different null rays 
by the null parameter 

u:=t-r*, (5.2) 

where r* is the "tortoise" coordinate 

r* :=r + 2mln(^- - l) . (5.3) 

These same null rays can also be labelled in terms of the proper time r of the vibrating 
trajectory. Given a vibrating trajectory we can obtain the relation u = u{U) we are 
searching for, if we identify U = t. By using the Schwarzschild metric it is straightforward 
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to see that 



u(U) = / dU 



l+fl-^j)" 1 ^) 



2 m 
r(U) 



1/2 



1/2 

r*(U) . (5.4) 



Now, given the relation u = u(U) one can calculate the Bogoliubov coefficients. As 
we are only interested in the presence of radiation with nearly Planckian spectrum, we can 
instead calculate directly the quantity [35, 36] 



d 2 U /df/V 1 d 2 u f du 



du 2 \du dU 2 \ dU 



(5.5) 

U(u) 



This function represents the peeling of geodesies in the geometry and is the magnitude 
relevant for encoding the radiative properties in situations that go beyond the standard 
scenarios usually analyzed in the context of Hawking radiation, for example, in the pres- 
ence of dynamical horizons or in the absence of a strict trapping horizon. As explained 
thoroughly in [35, 36], when this function is constant k(u) ~ k(u*) =: over a sufficiently 
large interval around a given u* , one can ascertain that during this same interval the system 
is producing a Hawking flux of particles with a temperature 

r=£. (5.6, 

Let us now calculate this quantity for the trajectories of interest. Instead of the 
variable r, let us use x := r — 2m. As we will always be in the regime x/2m C 1, we can 
simplify (5.4) which becomes u(U) ~ u a (C/), where 

u a (U) = 2m [ dU^ (l + V2 - 2mln JL (5.7) 



x 



so that 



2m J-L l + -2m-, (5.8) 

x \ 2mx z J x 



du a _ _ N , X \!/2 

dU ~ 

r \2 / • \ 2 1 I X _ II 

— - -2m sienfx) ( - I im±2 

, x\~ l XX x 

+ 2m 



x J \ x 2 

V 2mx 2 J 

Now, consider a trajectory x(r) such that it is in free-fall from x ou t to x- m where it bounces 
back symmetrically to return to x out (and so on and so forth). Except at the bounce itself 
the acceleration of the trajectory is 'x ~ — l/4m, i.e., it is equal to the surface gravity 
acceleration of the black hole. Then, if one plugs this value of the acceleration in (5.9), it 
is easy to realize that 

d 2 w a / dw a \ ~ 2 1 
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Figure 2. Depiction of the trajectory associated with a single vibration of a black star object. 
The trajectory is free-fall between x out — 0.0001 and Xf, = 0.00001, then it connects with an elastic 
bounce between x\> and x m = 0.000005. Finally, it goes back to x ou t in free-fall, before starting a 
new vibration cycle. We use 2m = 1 units. 

For instance, figure 2 shows a free-fall trajectory with a regularized bounce and figures 3 
and 4, the resulting behaviour for k(u). Therefore, a black star continuously vibrating in 
this way will continuously emit a Hawking-like flux of particles with a Planckian spectrum, 
only contaminated by a small non-Planckian emission associated with the bounces. 

Even if the vibration did not have these precise characteristics, one would still find a 
time interval within each cycle for which k(u) ~ constant and hence an overall approxi- 
mately constant k as long as the vibration cycle does not depart significantly from a free-fall 
trajectory. As seen from infinity the duration of a vibration half-cycle can be estimated to 
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Figure 3. Exact value of n(u) for the vibration of figure 2 in 2m = 1 units. 
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Figure 4. Exact value of k(u) for the vibration of figure 2 in 2m = 1 units (detail). The approximate 
k ~ (4m) _1 = 0.5 value is seen to be almost exact in the free-fall regions. 

be 

At v ~ 4m In— . (5.11) 

With a sufficiently large ratio x OVL t/xi n this duration can be made several times larger than 
Ky 1 = 4m, the inverse of the peak frequency of the Hawking spectrum. There would then 
be two options depending on the time resolution of the particle detector at infinity. With 
low temporal resolution we will see a Planckian spectrum with an additional peak at the 
vibration frequency, while with higher temporal resolution (but still low enough to see the 
Hawking spectrum), we will observe a Hawking spectrum with short periodic bursts of 
additional (non-thermal) radiation due to the bounces. 

In the case of the black stars of ref. [24], objects always on the verge of falling into 
their own gravitational radius, this would happen not only for the body as a whole, but 
also for each ball of matter contained within every internal radius. The pulsations of such 
a black star would thus resemble standard stellar pulsations in the sense of depending 
on the internal structure. Notice that, overall, the energy density profile of the black 
star would follow an almost l/(47rr 2 ) law. The closer the radius to the centre of the 
star, the bigger the surface gravity associated with this radius. Thus a vibration of an 
internal part of the black star will produce black body radiation with a larger temperature 
than the one associated with the surface gravity of the real external surface of the star. 
However, this radiation will be more red-shifted on its way out towards infinity. These 
two effects exactly compensate so that any radiation at infinity will be seen with the same 
temperature no matter at which precise radius it originated. The in principle non-coherent 
superposition of the radiation coming from vibrations at different radii would smooth out 
possible fluctuations on the value of k due to the actual vibration profiles. The final 
situation would be something practically indistinguishable from the perfectly continuous 
(non-fluctuating) emission deduced by Hawking. 
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In view of this discussion, it is not unreasonable that the generalized laws of black 
hole thermodynamics based on the geometrical notion of horizon area serve as a good 
approximation to the much "dirtier" physics of these material objects, which could have 
an area-entropy law and a Hawking emission at the appropriate temperature. 

Regarding the central theme of this paper, the point that we want to highlight is that 
the maximum frequency involved in the calculation of the particle production associated 
with these vibrations or bounces is not unreasonably large (remember that in this section 
we are assuming a relativistic dispersion relation) . The maximum frequency invoked in one 
vibration will be of the order 

uw ~ K H e KHA V (5.12) 

For a solar mass black hole and even for an oscillation with a semi-period as large as At v ~ 
100/kh, which corresponds to a ratio x out /x- m ~ 10 43 , this cut-off value u; max ~ 10 4 wp is 
perfectly compatible with the stringent bounds on Lorentz symmetry breaking imposed by 
current observations [37, 38]. To have a better idea of the behaviour of this magnitude, 
for an oscillation between x- m ~ 10 -2( Vh and x ou t — 10~ 10 rH, we have KuAt v ~ 23, so 
that uj max ~ 10~ 29 wp. When the system develops a new vibration, the frequencies involved 
in the process do not blue-shift further but are once again below the previous maximum 
frequency. There is then no trans-Planckian problem of any sort in this scenario. The 
maximum energy-scale for which we have to assume that a relativistic quantum field theory 
is reliable is perfectly reasonable. Therefore, within this scenario the prediction that there 
will be Hawking-like emission is absolutely reliable. Moreover, as we have shown, this 
radiation is Planckian in shape, but it is not thermal, in the sense that there aren't any 
hidden correlations to trace over. The partners of the Hawking particles are not trapped 
behind a horizon, so that they are also emitted towards infinity if only with a time delay. 
In the absence of a horizon, as well as of a singularity, the information-loss problem is 
completely absent in this scenario. 

6 Summary and conclusions 

The mechanism for thermal emission of a black hole found by Hawking contains a trans- 
Planckian problem that ultimately makes it extremely unreliable. Using this very problem 
as a guiding principle, we have presented and discussed three different scenarios, with 
three different Hawking-like emission mechanisms, in which the trans-Planckian problem 
does not show up. The trans-Planckian problem is caused by the combined effect of having 
a long-lived trapping horizon and a relativistic field theory at work up to huge energy 
scales. Thus, the three scenarios that we have analyzed are: 

1. Assume a long-lived horizon but effective modifications of the dispersion relations of 
subluminal type; 

2. Assume a long-lived horizon but effective modifications of the dispersion relations of 
superluminal type; 
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3. Keep a relativistic dispersion relation, but eliminate from the start the presence of 
strict long-lived trapping horizons. 

In the first scenario the existence of a long-lived trapping horizon at the classical level 
of general relativity (a unidirectional barrier for signalling) is preserved by the high-energy 
physics of the underlying theory. Therefore, this scenario suffers from the information- 
loss problem of the standard evaporation scenario and makes it impossible in practice to 
empirically explore the physics at work inside the horizon. 

In the second scenario the existence of a long-lived trapping horizon is just an effective 
description at low energies. The high-energy physics of the underlying theory allows for 
superluminal signalling and so the interior of the hole is exposed to the outside world. 
However, as we have shown, this scenario is not self-consistent in an astrophysical context. 
Under the reasonable assumption that superluminality would regularize the internal sin- 
gularity, configurations with a horizon become extremely unstable in astrophysical terms. 
Once formed, a Solar-mass black hole would tend to eliminate its horizon in a time scale of 
tens of microseconds. The apparently stationary astrophysical bodies standardly regarded 
as black holes should then be something different, without any horizons whatsoever. 

This leads us to the third scenario, in which we assume from the start that the astro- 
physical bodies usually called black holes are instead black stars, material objects hovering 
right outside but extremely close to their gravitational radius. We have shown that vi- 
brations of these objects can result in a Hawking- like emission. The maximum mode 
frequencies involved in this emission mechanism are very much under control and never 
reach large trans-Planckian values. Therefore, for the effective description of this emission 
mechanism, it is not necessary to consider any deviations from the relativistic dispersion 
relations. 

The emission of these objects will be Planckian although not strictly thermal. Since 
there aren't any horizons, there will be no hidden correlations and no information-loss 
problem: The Hawking partners will be emitted if only with a short time delay. On 
top of Hawking-like emission the spectrum will also show some short periodic pulsations 
associated with the bounces in each vibration cycle. This could eventually serve as an 
experimental signature of these objects. An effort to explore whether there exists a com- 
pelling dynamical mechanism behind the existence of these speculative black-star objects 
is on-going. 
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A Stability analysis 



In this appendix, we provide the details of the numerical analysis carried out to find 
the instabilities of a scalar field that obeys the modified two-dimensional Klein-Gordon 
equation 



(d t + d x v)(d t + vd x ) 



2 



-I* 



The dispersion relation associated with this equation is 

(w - kvf = k 2 (l + k 2 /kl) 



(A.l) 



(A.2) 



This dispersion relation is also found in BECs [22, 39]; in fact, the analysis below is similar 
to that of instabilities and quasi-normal modes in BECs [23, 40], except for using different 
boundary conditions. In particular, we will analyze flow velocities v with step-like profiles, 
i.e. such that they vanish in the exterior and have a constant and homogeneous value 
\v\ > 1 in the interior. By convention, the flow will be directed towards the left so that 
— 1 < v = tv < in the outer region x > while v = V[ < — 1 in the inner region x < 0. 
Furthermore, let us assume that the field is fully reflected at some point inside the hole 
located at a distance L from the horizon, i.e. at x = —L. 

The matching conditions at the point x = where the velocity profile is discontinuous 
can be obtained by requiring that the third spatial derivative of the field should also 
be discontinuous there and compensate for the jump in v in the generalized Klein-Gordon 
equation. This translates into the four conditions 



M = w\ = w 



o. 



[vd t 4> + v 2 cf)' + 4>'"/kl] 



0, 



(A.3) 



(A.4) 



where [/] := /(e) - /(-e). 

For a given frequency u), the general solution can be written as 

' 4 

3=1 
8 

J^Aje^-^ (x>0), 
U= 5 

where {kj} are the roots of the corresponding dispersion equations (four roots for each 
homogeneous region), and the constants Aj have to be such that the matching conditions 
(A.3) are satisfied. We can write down these conditions in matrix form Aij(ui)Aj = 0, 
where the transpose of A is 



/ 


1 


fcl 


1 


VlU - 


vfh + kf/k 2 \ 




1 






VlU - 


vfk 2 + kl/kl 




1 


h 


^3 


VlU - 


vfk 3 + kl/kl 




1 


k,{ 


kf 


VlUJ - 


v^k^. -\- k^/kp 




-1 


-k 5 




— {v t oj 


- v 2 k 5 + kl/k 2 ) 




-1 


-ke 


-k 2 


— {v t oj 


- v 2 k 6 + kl/k 2 ) 




-1 


-k 7 




— {v t oj 


- v 2 k 7 + k^/kl) 


V 


-1 


-kg 


-k 2 
K 8 


— (v r u) 


-v 2 k s + kl/k 2 ,)) 



(A.5) 
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Furthermore, to find the unstable modes these conditions have to be complemented with 
conditions at the boundaries of the system, namely at x — > oo and x = —L. We will deal 
only with instabilities so we will be concerned with modes with Im(uj) > 0. Furthermore, 
for each mode with frequency u, there is another mode with frequency —a;*. So we will 
consider only Re(w) > without loss of generality. 

For complex oj the four roots of the dispersion relation in the right region are complex: 
Two of them with positive imaginary part (say k§ and ke) and two of them with negative 
imaginary part (say k 7 and ks) [22]. Since the modes must be normalizable, we must ensure 
that 



A 7 = As. 







(A.6) 



so that the solution does not grow exponential at x — > oo. No additional requirements (e.g. 
outgoing conditions) are imposed at infinity so the boundary conditions we use match the 
analysis in [16, 22, 41]. 

Let us now analyze the reflection at x = — L. In order to find appropriate boundary 
conditions, we have to notice that in the presence of dispersion, the current 



1 



n, I ) 



(A.7) 



is locally conserved, so that it should vanish at x = —L. One can choose different reflect- 
ing boundary conditions (generalization of the Dirichlet, Neumann and Robin boundary 
conditions). Qualitatively our results do not depend on this choice. For concreteness we 
have used generalized Dirichlet conditions: 



4>{-L) = 0, cj>'{-L) = 0, 
which amount to the following conditions on the modes: 



(Ai 



0. 



e — ikjL 



0. 



(A.9) 



These four conditions (A.6) and (A.9) can be imposed by adding the following extra 
rows to the 4x8 matrix A(w) : 



/0 0001 o\ 

1 

e -ifciL e -ifc2-L e ~ik'iL e ~ikiL q q q q 

\kie~ iklL k 2 e~ ik2L k 3 e~ ik3L k A e~ ikiL 0/ 



(A.10) 



The extended 8x8 matrix A(uj) thus obtained is such that the possible unstable modes 
must satisfy ^ • Aij(oj)Aj = 0. This set of equations will have a solution for every value of 
oj at which its determinant vanishes. 

The numerical method then follows these steps: 
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1. Given the values of the flow velocity in the left v\ and the right v r , for each frequency 
a; in a grid covering the first quadrant of the complex plane, we calculate its associated 
/c-roots [by solving the dispersion relation (A. 2)] in both regions. 

2. Given the value of the size L of the interior of the black hole, for each value of oo, 
we evaluate the determinant det[A(w)] . The values of oj for which this determinant 
vanishes correspond to instabilities of the system. 

3. Instead of drawing contour plots of det[A(u;)], in order to increase contrast, we plot 
the following function: 

F(u) = - log[/(u;)]/max{- log[/(w)]}, (A.ll) 

where 

f(uj) = g(u)/max[g(uj)}, 

g{uj) = det[A(w)]/([det A(w)]), (A.12) 

and (h) stands for a coarse graining of the function h that removes the peaks, hence 
leaving only the background behaviour of h. The maxima of this function F(oj) 
correspond to vanishing values of det(A) and hence to instabilities of the system, as 
in figure 1 above. 
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